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I. INTRODUCTION 



The interior of neutron stars is, to a very good approximation, formed by pure neutron 
matter [1, 2]. At the very initial stages after their formation, these objects are very hot, 
with temperatures as high as T ~ 40 MeV [3]. The Equation of State (EoS) of pure neutron 
matter in a wide range of densities and temperatures is therefore a crucial ingredient to de- 
scribe the structure and the evolution of neutron stars. The evaluation of both the neutron 
matter and the symmetric nuclear matter EoS starting from realistic models of the nucleon- 
nucleon (NN) interaction is still a major challenge in nuclear physics. The short-range and 
tensor components of realistic NN forces induce correlations which substantially modify the 
many-nucleon wave function as compared to the free Fermi gas (FFG) Slater determinant. 
This is particularly important for symmetric matter, where the ^Si-'^'Di channel plays a 
pivotal role. In neutron matter, Pauli effects block this tensor channel, but short-range 
correlations still need to be accounted for appropriately. Several theoretical approaches 
have been developed over the years to treat these correlations in zero temperature neu- 
tron matter: variational techniques within correlated basis functions [4-6]; Auxiliary Field 
[7] or Quantum Monte Carlo [8] calculations with simplified interactions and the popular 
Brueckner-Bethe-Goldstone hole-line expansion [9] in its lowest order form, the so-called 
Brueckncr-Hartree-Fock (BHF) approximation [10]. At finite temperatures, fewer efforts 
have been focused in this direction: the well-known variational calculation of Friedman and 
Pandharipande [11] and recent similar calculations [12], as well as BHF extensions at finite 
temperature [13, 14]. The latter approximation takes into account particle-particle corre- 
lations by solving the Bethe-Goldstone equation, which leads to the so-called G-matrix. 
Nevertheless, a minimal consistent treatment of correlations in nuclear systems requires the 
inclusion not only of particle-particle (pp) intermediate states, but also of the hole-hole (hh) 
ones. The propagation of particles and holes can be treated in the same footing by means 
of the Self-Consistent Green's Function (SCGF) approach [15]. 

The SCGF approach gives direct access to the single-particle spectral function and there- 
fore to all the single-particle properties of the system. A great progress in the application 
of the SCGF method to nuclear matter has been achieved in recent years, both at zero [16] 
and finite temperatures [17-21]. The solution of the SCGF equations is a rather demanding 
numerical problem due to the complete treatment of off-shell energy dependences. As a 
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consequence, the SCGF method has been apphed to few general, extensive analysis of dense 
nuclear systems. The studies at zero temperature have been mainly oriented to provide the 
appropriate theoretical support for the interpretation of (e, e'p) experiments, while those at 
finite temperature focus on a correlated description of matter to be used in the studies of 
heavy ion collisions dynamics or in astrophysical environments. In particular, the effects of 
temperature might affect substantially different astrophysical observables. As an example, 
the cooling curve of a neutron star depends on the interior temperatures and the possible 
transition to a superfluid regime [22]. Also, the gravitational wave signature of the super- 
nova explosion might be sensitive to the EoS and might even be able to distinguish thermal 
effects [23]. 

In this line, we want to study the microscopic and thcrmodynamical properties of hot 
pure neutron matter within the SCGF framework. The SCGF method, as formulated here, 
cannot be used below the critical temperature of the pairing transition [24, 25] and therefore 
all our results only apply for the normal phase. Although this is not the first time that the 
SCGF approach is used to study pure neutron matter [17, 26], it is, up to our knowledge, the 
first time that a systematic study of the microscopic and thcrmodynamical properties of pure 
neutron matter at finite temperature is performed within the SCGF approach. Moreover, we 
shall perform our calculations with two different realistic nucleon-nucleon interactions, the 
meson-exchange CD-BONN potential [27] and the local Argonne V18 [28]. Together with 
the comparison to other many-body approaches, this can be used to highlight the model 
dependence in hot neutron matter calculations. 

Lately, the problem of neutron matter has also been growing in interest due to its connec- 
tion with the experimental studies of ultracold fermionic systems [8, 29]. Dilute strongly- 
interacting fermionic systems with large scattering lengths (such as neutron matter, with a 
scattering length a — —18 fm to be compared to kp — 1.68 fm~^ for p — 0.16 fm~^) lie in 
the so-called unitary regime. As a consequence of the lack of any characteristic energy scale, 
these systems show a universal behavior in their zero- and finite-temperature dynamics, with 
scahngs that are related to the non-interacting case [30]. We shall not treat this particular 
problem here, but one should mention that the SCGF method is able to tackle the unitary 
regime above the pairing phase transition [31]. When properly complemented with pairing 
effects [24, 25], this method should also be able to properly describe the unitary regime. 

After a brief description of the SCGF formalism in Section II, we discuss in Section III 
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our results for the microscopic properties of hot pure neutron matter. Section IV is devoted 
to the analysis of the thermodynamical properties and the comparison of our results with 
those obtained within other approaches. Finally, a brief summary and our main conclusions 
are presented in Section V. 



II. SELF-CONSISTENT GREEN'S FUNCTIONS METHOD AT FINITE TEM- 
PERATURE 

A crucial step in the microscopic description of nuclear many-body systems is the de- 
termination of the effective in-medium nucleon-nucleon (NN) interaction. The ladder ap- 
proximation to the in-medium T-matrix is well suited for strongly interacting low density 
systems [32] and has the following structure: 

(kik2|r(Q+)|k3k4) = (kik2|y|k3k4) 

+ I |^^(kik2|l-|k5k6)^°,(A;5,A;6;0+)(k5k6|r(0+)|k3k4) , (1) 

where Qjj is associated to the propagation of two dressed but non-interacting single-particle 
lines: 

J-oo ~ ^ ~ ^ 

with f{u!) — ^e^i'^~^^^ + l] ^ the Fermi-Dirac distribution and A{k, cu) the single-particle 
spectral function. The notation Q± stands for Q ± irj, with rj infinitesimally small. Q^j can 
be interpreted as a Pauli blocking factor at finite temperature, analogous to the one that 
appears in zero temperature BHF calculations [15]. In contrast to BHF, however, the zero 
temperature version of the SCGF formalism accounts for the intermediate propagation of 
both pp and hh states. 

The interaction of a nucleon with the remaining nucleons in the medium is described 
within the Green's functions formalism in terms of the self-energy [33] . Its imaginary part 
is related to the in-medium T-matrix: 

|^<kk'|Imr(a; + u;;)|kk'>^(A;,u;')[/(^')+K^ + ^')], (3) 

where a Bose- Einstein factor, = [e^/^C^-^M) _ ^ ^ appears due to the symmetric treat- 
ment of pp and hh states. The real part of the self-energy is determined from its imaginary 
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part by a dispersion relation: 

ReE{k, u) = EHF{k) - V / , (4) 

J TT 00 — OU 

except for the energy-independent Hartree-Fock contribution: 

^HFik) = / ^ (kk'l V|kk') nik') , (5) 

where the momentum distribution includes the effects of correlations via A{k,u!): 

n{k)^i. 1^ ^A{k,u;)f{u;). (6) 

1/ — 2 accounts for the spin degeneracy of neutron matter. Finally, one can make use of 
Dyson's equation to close this set of equations by determining the single-particle spectral 
function from the real and imaginary parts of the self-energy: 

Aik,uj]^ 5 -. (7) 

[^-^-ReE(^,a;)]' + [ImE(^,^)f 

The previous equations are derived within the grand-canonical picture, where the two 
external, fixed variables are the temperature, T = ^, and the chemical potential, fi. For 

dense matter studies, it is more convenient to fix the density p and therefore we supplement 
the previous set of equations with the normalization condition: 

f d^k p duj 

which determines a "microscopic" chemical potential, /2. In a thermodynamically consistent 

approximation (such as the ladder approximation), jl should coincide with the macroscopic 
chemical potential, fj,, obtained from the bulk properties by taking the derivative of the free 
energy density, F: 

dF 



dp 



(9) 

T 



Thermodynamically non-consistent many-body approximations, such as BHF, lead to p, ^ p 
[34]. 

Equations (1-8) form a closed self-consistent set of equations in terms of the in- medium 
interaction, the self-energy and the single-particle spectral function that can be solved it- 
eratively. The numerical details associated to the solution of these equations are rather 
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involved and we refer the reader to Refs. [19, 35, 36] for further details. It is important 
to note that the numerical solution of the SCGF method, when available (see following 
paragraph), accounts for the full ladder approximation. 

The bosonic factor appearing in Eq. (3) presents a pole for Q = 2/i, which is generally 
cancelled by an associated zero in ImT(0 = 2/x). However, below a certain critical tempera- 
ture, Tc, the state with center of mass momentum P — and energy fl — 2fj, does not cancel 
the bosonic factor and an instability occurs, reminiscent of the formation of Cooper pairs. 
This signals the onset of superfluidity, according to the so-called Thouless criterion [37, 38] , 
and imposes a limit to the lowest temperatures we can achieve within our numerical calcu- 
lations. All the results presented in the following are obtained for T > thus neglecting 
the effect of pairing correlations but guaranteeing the convergence of the approach. 

So far, we have discussed the determination of the microscopic properties of the system. 
The Green's function formalism can also be used to obtain the bulk properties of neutron 
matter. For the case of two-body interactions, one can evaluate the energy per particle by 
means of the Galitskii-Migdal-Koltun (GMK) sum rule [39, 40]: 

7^::<7—+uj\A{k,uj), (10) 



A p J (27r)3 271 2 {2m 

from the spectral function evaluated in the SCGF approach. To obtain the free energy, 

F = E — TS, and have a complete thermo dynamical description of the system, one still 

needs to compute the entropy within a correlated approximation. This can be obtained by 

using the Luttinger-Ward (LW) formahsm [21, 41, 42]. Within this approach, an expression 

for the grand-canonical potential in terms of dressed single-particle propagators can be 

obtained by means of a Legendre transformation. The entropy can then be computed from 

the derivative S = — |p |^, which gives a closed expression in two terms, S = S^^ + S'. The 

first one corresponds to the dynamical quasi-particle (DQ) entropy density: 

Sdq V f d^k duj 

— ' ' — a{u!)B{k,u!) , (11) 



A pj {2nfJ_^2Ti 
given by the convolution of a statistical factor, (t{uj) — —f{oj) In /(a;) — [1 — /(<^)] In [1 — /(<^)]) 
and a spectral function, B{k,uj): 



B{k,uj) 



_ dJ:{k,oj) 



A{k, uj) - 2^Mi^ImE(;t, uj) , (12) 



which can be computed from the single-particle quantities obtained in the SCGF approach. 
This B— spectral function accounts for the effect of the dynamical {i.e. interaction- induced) 
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correlations that fragment the quasi-particle peak [21]. In this paper, we will consider 
that the second term, S', is negligible due to constraints in phase space for relatively low 
temperature [42]. This approach leads to thermodynamical consistent results for neutron 
matter as well as for symmetric nuclear matter [21, 36]. 

In order to assess the dependence of our results on the many-body approximation em- 
ployed in the description of neutron matter, we shall compare the SCGF calculations to a 
finite temperature generalization of the BHF method. A real finite temperature extension of 
the BHF approach is given by the Bloch-de Dominicis theory [43, 44] but, instead of using 
the latter, our calculations will rely on an often used simpler generalization [14, 45]. This 
extension can be obtained from the SCGF equations by assuming that the spectral function 
has no width and full strength concentrated at the BHF quasi-particle energy: 

Aik, Lu) = {27i)6 [a; - eBHF{k)] . (13) 

In addition, one eliminates the bosonic factor of Eq. (1) and modifies the in-medium two- 
body propagator to include only intermediate particle-particle propagation: 

\l+ - SBHFKk) - SBHFKk') 

The set of equations thus obtained mimics the zero temperature BHF formalism with the 
replacement of the step-function momentum distributions at T = by Fermi-Dirac distri- 
butions at T 7^ 0. This guarantees that in the T ^ hmit the results will coincide with 
BHF. One can proof that this extension coincides with the Bloch-de Dominicis results at 
low temperatures [44]. 

Few other approaches exist that can be used to study neutron matter at finite tempera- 
tures starting from reahstic NN potentials. The benchmark variational calculations of Fried- 
man and Pandharipande (hereafter FP) [11] relied on a frozen correlation approximation, 
i.e. using as a starting point the Jastrow-like correlation functions obtained at zero temper- 
atures. This is of course an additional approximation, possibly only suitable for low enough 
temperatures and large densities, where matter can be consider degenerate. The variational 
approach has only recently been extended to finite temperatures to include appropriately 
thermal correlations [46] . Alternatively, the case of low densities and high temperatures can 
be studied by means of the model- independent virial expansion [47, 48] . In this approxima- 
tion, the thermodynamical properties are expanded in terms of the fugacity, z — e^^. The 



first term in this expansion leads to the thermodynamics of a classical free gas, while the first 
order correction is given in terms of a virial coefficient that can be computed from the exper- 
imental NN interaction phase-shifts in free space. Since neutron matter is not expected to 
clusterize at low densities, this approximation will hold for extremely dilute and hot matter. 
Recently, another method has been proposed to study neutron matter at nonzero temper- 
atures by making use of renormalized low momentum two- and three-nucleon interactions 
whose short-range components have been properly eliminated [49]. The thermal properties 
of neutron matter have been computed up to second order in finite temperature many-body 
perturbation theory, including contributions from normal and anomalous diagrams. 

III. MICROSCOPIC PROPERTIES OF NEUTRON MATTER 

In this Section we will discuss the microscopic single-particle properties of neutron matter 
as obtained from the SCGF approach. To address the model dependence of our calculations, 
we will show results using two different realistic nucleon-nucleon interactions, namely, the 
meson-exchange CD-BONN [27] and the local Argonne V18 potentials [28]. Partial waves 
up to J = 8 have been considered, with the Born approximation for J > 5 in both SCGF 
and BHF calculations. 

We start by showing in Fig. 1 the density and temperature dependence of the neutron 
spectral function in dense neutron matter. Due to the similarity of the results for the two 
interactions, we will only consider the results obtained with the Argonne V18 potential. 
The spectral function for densities ranging from p — 0.04 fm~^ to p = 0.32 fm~^ at a fixed 
temperature of T = 5 MeV is shown in the left panels for three momenta: = (top 
panel), k = kp (middle panel), and k = 2kp (bottom panel), kp corresponds to the Fermi 
momentum associated to each density. The right panels show the results for a fixed density, 
p — 0.16 fm~^, and temperatures from T = 5 to 20 MeV. The qualitative features of these 
figures are already well-known (see e.g. Ref. [35]). There is an important quasi-particle 
peak, which contains roughly 70 — 80% of the total strength for all momenta. The position 
of this peak changes with momenta and it is described by the self-consistent equation: 

Sqpik) = + Re E[A;, eqp{k)\ , (15) 

which defines the quasi-particle spectrum for neutrons in the medium. With increasing 
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density, the quasi-particle peak at zero momentum shifts to lower energies with respect to 
the chemical potential. It turns out that neutrons at low momenta are more bound at higher 
densities. The situation is the opposite for high momenta {k ~ 2^^?), where the peak shifts 
to higher energies when density increases. At the Fermi surface, k = kp, the quasi-particle 
peak is approximately centered around u ^ /j, and its width decreases as p increases. At 
zero temperature and in the absence of pairing correlations, the spectral function would 
actually have a delta-like quasi-particle peak. The effect of density is particularly large in 
the low- and high-energy tails of the spectral function. For both large removal {to « ji) 
and large addition {uj » fx) energies, the strength increases with density. These off-shell 
components of the spectral function are populated mainly due to the action of the short- 
range core of the nuclear interaction and therefore it is reasonable that they increase when 
the mean separation between neutrons decreases. In other words, the high-energy strength 
of the spectral function is a good measure of the correlations induced by density effects. 

The influence of temperature in the spectral function is less pronounced. Both the posi- 
tion of the quasi-particle peak and the strength at low and high energies are almost unaf- 
fected by changes in temperature. The only region that is slightly modified by temperature 
corresponds to the range of energies ~ yU, which is particularly sensitive to variations in 
phase space [50] . It seems fair to say that the structure of the spectral function is mainly 
determined by the in-medium renormalization associated to the density, while temperature 
effects play a minor role. This is no longer the case close to and below Tc, where a rela- 
tively small decrease in temperature can lead to the appearance of superfluidity and thus to 
an important change in the properties of the spectral function. In particular, the onset of 
pairing results into a double quasi-particle peak structure close to the Fermi surface [24, 25]. 

To learn more about the effect of hh propagation on the microscopic properties of neutron 
matter, one can compare the quasi-particle peak described by Eq. (15) with the single- 
particle spectrum obtained within the BHF approach. The first includes both pp and hh 
effects, while the second only accounts for pp states. In Fig. 2 we compare the real part 
of the on-shell self-energy, ReT^lk, eqp{k)], for both approaches at densities p = 0.08, 0.16 
and 0.24 fm~^ (left, central and right panel, respectively) and temperatures T = 5 MeV 
(solid lines) and T = 20 MeV (dashed lines). The results displayed in this Figure have been 
obtained with Argonne V18, but similar conclusions are reached with CDBONN. For all 
cases, the SCGF spectra are more repulsive than the BHF ones at all momenta. The effect 



9 



of hh propagation in the on-shell self-energy is therefore of a repulsive nature. This effect is 
larger at low momenta, in accordance with the idea that the dressing induced by hh states 
is irrelevant for high-momentum, particle states. The repulsive effect of SCGF with respect 
to BHF increases with density and the differences can be as large as 25 MeV for /c = at 
p = 0.24 fm-^ 

The temperature behavior of the quasi-particle spectra shows some interesting features. 
On the one hand, the BHF single-particle spectrum becomes more repulsive with increasing 
temperature at all momenta. This is usually attributed to the presence of thermal Fermi- 
Dirac factors in the self-energy. The repulsive high relative momentum components of the 
interaction are not accessible at zero temperature and they only become available once the 
thermal distribution populates high momentum single-particle states. The overall effect is 
then repulsive. The same reasoning applies to particle states in the SCGF case, which also 
become more repulsive with increasing temperature. Hole states, on the other hand, become 
more attractive with increasing temperature. Presumably, this behavior can be attributed 
to the fact that hole states are renormahzed in the SCGF, which results into a quenching 
of the attractive long-range components of the NN interaction in the zero temperature case. 
The inclusion of thermal effects leads to a somewhat weaker renormalization that increases 
the attractive component of the spectrum for k < kp- A similar effect has been observed in 
extended BHF (where the repulsive contribution of holes is taken into account by the M2 
rearrangement term in the self-energy [51]) as well as in SCGF calculations of symmetric 
nuclear matter [36]. 

Among the one-body properties of interest for correlated many-body systems, the mo- 
mentum distribution of Eq. (6) is particularly sensitive to dynamical corrections. At zero 
temperature, for instance, the momentum distribution of the FFG is just a step-function, 
with complete population below kp and empty states above. In contrast, the correlated 
momentum distribution at T = displays a substantial depletion of hole states and a non- 
zero population of high momentum states. Unfortunately, the FFG at finite temperature 
also shows these features, since all the states become partially populated due to the thermal 
distribution of states. As a consequence, the correlated n{k) at finite temperature will have 
both thermal and dynamical components. To appropriately disentangle these two com- 
ponents, extensive studies of the temperature and density dependence of the momentum 
distribution are needed. 
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This analysis is presented in Fig. 3, where the density (top panels) and temperature 
(bottom panels) dependence of the momentum distribution is shown for both CD-BONN 
(left panels) and Argonne V18 (right panels) potentials. Interesting analogies between the 
density and temperature dependence are observed: decreasing temperature has a somewhat 
similar effect to increasing density. This is in stark contrast to the effect of density and 
temperature on the spectral function, which, as already commented, are rather different. In 
the case of the momentum distribution, these dependences can be interpreted in terms of 
degeneracy arguments: for both the low temperature and the high density case, the system 
is reaching a degenerate limit, where thermal effects are unimportant and the depletion is 
essentially governed by dynamical effects. This will be the range which is interesting for 
understanding the influence of short-range correlations on the system. The opposite limit 
(low densities, high temperatures) leads to a momentum distribution which is controlled by 
thermal effects. 

Some particular details, however, differ depending on how the degenerate limit is ap- 
proached. On the one hand, fixing the density and progressively decreasing the tempera- 
ture, leads to a monotonous increase (decrease) of n{k) below (above) the Fermi surface. In 
particular, n{k = 0) saturates to a value different from 1 when T ^ 0. For p = 0.16 fm~^ 
at the lowest temperature available {T — A MeV), one finds n(0) = 0.974 for CDBONN and 
n{0) — 0.959 for Argonne V18. The differences in the short-range components of the two 
interactions explain the discrepancies in n{0): Argonne V18 has a harder short-range core 
compared to CD-BONN and thus leads, in general, to lower occupations ioi k < kp Sit high 
densities. On the other hand, fixing the temperature and increasing the density, one finds a 
different scenario, where n(0) is no longer a monotonous function due to the competition of 
thermal and dynamical effects. 

This behavior is observed in detail in Fig. 4, where the occupation of the lowest mo- 
mentum state, n(0), is shown as a function of density for several temperatures. The density 
dependence of n{0) indeed has features which can be attributed to both thermal and dynam- 
ical effects. For all temperatures, there is a steep decrease of n(0) when p ^ 0. The FFG 
n(0), shown in double-dotted dashed line, has a similar behavior, which can be explained 
in terms of the system approaching the classical limit (/i — > oo). In the non-interacting 
case, dynamical correlations are absent and therefore thermal effects are responsible for the 
strong decrease of n{0) at low densities. The analogous behavior in the correlated n{0) is 



11 



basically driven by thermal correlations. The high density behavior of n(0), on the other 
hand, is totally different from the FFG. While the latter always equals 1 above p ~ 0.08 
fm~^ at T = 5 MeV, the correlated n(0) at this temperature tends to have values which are 
about 10% lower. One actually observes a decrease in n(0) as density increases in this low 
temperature range. This dependence can be understood in terms of dynamical correlations: 
an increase in density results into a decrease of the mean distance between particles and, as 
a consequence, the importance of short-range effects is incremented at higher densities. As a 
consequence, the depletion increases with density, as observed. Again, this effect depends on 
the particular short-range structure of the NN force, which explains the differences observed 
between the left and right panels. Finally, let us note once again that the temperature de- 
pendence of n(0) is monotonous: large temperatures lead to low n(0)'s at all densities. The 
changes induced by temperature on n(0) are however density dependent and, as expected 
from degeneracy arguments, they are almost negligible at high densities. 

IV. THERMODYNAMICAL PROPERTIES OF NEUTRON MATTER 

The SCGF approach, complemented with the Luttinger-Ward formalism, can be used to 
obtain the thermodynamical properties of neutron matter including the effect of correlations. 
In this Section we shall analyze these properties and compare the SCGF results with those of 
other approaches, such as the variational calculation of FP, the finite temperature extension 
of BHF and the virial expansion. 

The energy per particle, obtained from the GMK sum-rule of Eq. (10), is shown in 
Fig. 5 as a function of density for two temperatures, T = 10 and T = 20 MeV. CDBONN 
(Argonne V18) results are displayed in the left (right) panel. The SCGF results (circles) 
are compared with those obtained with the finite temperature generalization of the BHF 
approach (triangles) , and those of the variational calculation of FP (crosses) . Note that the 
results for the energy per particle are not quoted in the original publication and these have 
been obtained from the free energy and the entropy. At low densities, we also compare our 
results with the model-independent virial approximation for fugacities up to z = 0.5 [48]. 
These correspond to densities p = 0.0035,0.0098,0.0181 and 0.0279 fm~^ at temperatures 
T = 5, 10, 15 and 20 MeV, respectively 

Comparing the SCGF and BHF approaches for a single NN interaction, one finds that 
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the inclusion of hh correlations leads to a more repulsive energy per particle for almost 
all densities. As expected from phase space considerations, this repulsive effect is more 
important at higher densities. Moreover, the repulsion induced by hh propagation is more 
important for Argonne V18 than for CDBONN. As mentioned previously, the Argonne V18 
interaction has a strong short-range core and therefore the hh renormalization on top of the 
pp propagation will still have an important effect. In particular, at a temperature of T = 10 
MeV, the inclusion of hh propagation leads to a 1.6 MeV (4.0 MeV) increase of the energy 
per particle at p = 0.16 fm~^ (0.32 fm~^). In contrast, the weaker short-range structure of 
CDBONN is already well treated with pp correlations and the inclusion of the hh component 
has a smaller effect, of only 0.6 MeV (1.5 MeV). These results are in agreement with the zero 
temperature calculations of the Ghent group, which showed almost no difference between 
SCGF and BHF at p = 0.16 fm~^ for the Reid93 interaction [26]. However, these findings 
seem to disagree with those of the Krakow group [52], which suggest differences between 
the SCGF bulk energies and continuous choice BHF calculations of about 5 MeV in the 
same conditions. Note, however, that those results were obtained with a simpler separable 
NN interaction and that different numerical procedures were used in the solution of the 
SCGF equations. The recent calculation of Ref. [49] leads to more repulsive results than 
ours at low densities, even when three-body effects are not considered. This is curious since, 
by construction, Viowk does not include short-range cores and thus one would have naively 
expected theirs EoS to be softer than that obtained by renormalizing interactions with hard 
cores. 

The differences in energy between the many-body approaches for a given potential that 
we have just discussed are a consequence of differences in the treatment of dynamical and 
thermal correlations. In contrast, the discrepancies within the same many-body approach 
for two NN interactions are a reflection of the different structure of the two potentials and, 
in particular, of their short-range behavior. In general, the results with Argonne V18 for 
both the SCGF and the BHF approaches are more repulsive than those of CDBONN. In the 
SCGF approach at T = 10 MeV and 0.16 fm-^ for instance, E/A = 18.0 MeV for Argonne 
V18, while E/A = 16.9 MeV for CDBONN. This discrepancy increases with density and, 
for p = 0.32 fm~'^, it becomes as large as 6.4 MeV. In contrast, the differences in energy 
per particle between the two potentials for the BHF approximation are rather small. For 
T = 10 MeV and at p = 0.16 fm^^ they are less than 0.5 MeV, while at p = 0.32 fm-=^ 
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they are just about ~ 3 MeV. This indicates that the inclusion of hh correlations in the 
energy per particle increases the dependence of the results on the short-range structure of 
the potential. Let us also note that the discrepancy in the energy per particle due to the 
use of different NN potentials is somewhat larger than that associated to the use of different 
many-body approaches, particularly in the high density regime. 

At low densities, short-range effects are weakened and the SCGF data agrees with the 
virial expansion independently of the NN interaction. This is particularly well observed in 
the inset of Fig. 5. This agreement provides for the first time, to the best of our knowledge, 
a model-independent verification of the numerics of the SCGF approach. Let us also note 
the relatively large differences between the FP and the SCGF results below 0.08 fm^'^. In 
addition, we would like to stress that the various approximations reach the correct classical 
limit, E/A —>■ 3T/2 for p ^ 0, but the way this limit is reached depends on the approach 
under consideration. In all cases, the energy per particle shows a well defined minimum. 
This is a consequence of the competition between thermal effects, which are dominant at 
low densities and tend to make the energy more repulsive, and interaction effects, which are 
attractive and important at intermediate densities. 

Remarkably, our SCGF results for Argonne V18 agree well with those of FP at high den- 
sities. Both calculations are based on local NN potentials, but the Urbana V14 interaction 
of Ref. [11] includes a density-dependent quenching of the two-pion exchange contribution 
to account for the repulsive effect of a three-body force in a phenomenological way. Naively, 
one would have expected the inclusion of such a contribution to yield more repulsive results 
than ours, especially at high densities. Note that, if the contribution of the three-body ef- 
fects was indeed negligible, the observed agreement could be a signature of an unprecedented 
agreement between the variational and SCGF approaches in a wide range of temperatures 
and densities. To clarify this issue, it would be interesting to compare our SCGF results with 
finite temperature variational calculations with the Argonne VIS interaction. Alternatively, 
we have performed some preliminary SCGF calculations with the Urbana V14 interaction 
(together with the density-dependent quenching) . Our preliminary calculations indicate that 
the energy per particle is ~ 3 MeV more repulsive than the SCGF energy with Argonne VIS 
at p = 0.16 fm~^ and T = 10 MeV. The discrepancy increases to ~ 10 MeV at 0.32 fm~^. 
All in all, this seems to indicate that the agreement between the FP calculations with the 
Urbana V14 force and our SCGF with Argonne VIS is a coincidence, which might have been 
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caused by some sort of cancellation between the differences induced by the two underlying 
interactions and those associated to the different many-body approaches. 

The discrepancies are substantially smaller in the case of the entropy per particle, shown 
as a function of density for two temperatures, T = 10 and T = 20 MeV, in Fig. 6. In 
particular, the changes arising from the use of different potentials (left and right panels) 
are smaller than those due to the use of different many-body methods. At T = 10 MeV, 
the different approximations (DQ entropy from SCGF results, BHF entropy, FP, FFG) 
are quite consistent with each other. The deviations above 0.16 fm~^ are at most of 0.15 
Boltzmann units, which would have a maximum impact on the free energy per particle of 
T X SS/A ~ 1.5 McV. At higher temperatures (T = 20 MeV), the differences between 
approaches are somewhat larger, of at most 0.25 units. 

All in all, these results support the idea that the entropy is mostly determined by thermal 
correlations and rather unaffected by dynamical correlations. This is confirmed by the 
extremely narrow quasi-particle peak of the B spectral function. The many-body effects 
that fragment the quasi-particle peak, which are extremely important in the calculation 
of the energy, are almost negligible for the entropy [21, 36]. This explains partially the 
good agreement between SCGF and BHF entropies. The latter are obtained by using the 
quasi-particle approximation to the entropy: 



Although, as observed in Fig. 2, the quasi-particle energies of the two approaches are quite 
different, the change in chemical potential between BHF and SCGF shifts the entropy to 
values very close to those of Sdq- The similarity between both entropies had already been 
observed for symmetric matter [21]. 

Compared to the FP entropy, we find that both SCGF and BHF predict shghtly larger 
entropies at large densities for both temperatures and interactions. A similar effect was 
observed in Ref. [49] and attributed to an anomalously low effective mass in variational cal- 
culations. The restriction to a quadratic spectrum in variational approaches is a limitation, 
especially in view of the clearly non-quadratic momentum dependences of the BHF and 
SCGF quasi-particle spectra (see Fig. 2). The entropies in the latter approaches go beyond 
such approximation and, in the SCGF case, they even go beyond the assumption of a single 
quasi-particle peak. In any case, close to the degenerate limit, all the calculated entropies 




(16) 
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have a Fermi-liquid-like behavior: 

J-^sT, (17) 

(see Fig. 10), where the parameter = -^p^ is given in terms of the effective mass m*, 
calculated at the Fermi surface at zero-temperature. The discrepancy between the S^'^ and 
the FP entropies at large densities (close to the degenerate hmit) suggests that the effective 
mass in the DQ entropy is larger than that of the variational entropy. Indeed, the m* in 
the DQ density of states is the product of the ml and the m* effective masses [53]. The 
latter is associated to the energy dependence of the self-energy and is believed to be absent 
in the variational approach. Since is strongly peaked around the Fermi-surface, it leads 
to larger values of the total effective mass and therefore increases the DQ entropy at large 
densities with respect to the variational one. 

In this direction, it is important to note that, in all cases, the entropy at high densities 
is smaller than that predicted by the FFG. This is in accordance with the idea that, in this 
regime, the entropy is dominated by the effective mass, which is always smaller than 1 and 
therefore leads to lower entropies. Finally, let us stress that, as expected, the interaction 
has little influence in the entropy near the classical regime. At low densities, all the approx- 
imations to the entropy converge to similar values and no differences are observed between 
the FFG and the virial entropies (see inset of Fig. 6) . 

The free-energy obtained from the GMK sum-rule complemented with the dynamical 
quasi-particle entropy is shown in Fig. 7 as a function of density for several temperatures. Let 
us first note that the calculations yield well-behaved results in a large range of densities and 
temperatures for both CDBONN (left panel) and Argonne V18 (right panel). In particular, 
the low density high-temperature regime agrees well with the virial results. This agreement 
is directly related to the similarity of the DQ and the virial entropies since, in this regime, the 
entropy overcomes the energy contribution in F/A. Comparing the two panels, one observes 
that for densities higher than 0.08 fm~^ the Argonne results are more repulsive than the 
CDBONN ones. In addition, the Argonne SCGF and the FP results are quite close to each 
other for all densities, with differences (mostly coming from the entropy) smaller than 3 
MeV for the highest density considered here. In general, one can say that, for low densities, 
F/A is well determined and all the approaches agree well with each other independently of 
the potential. Above ~ 0.08 fm~^, however, differences appear due to the sensitivity of the 
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many-body approach to the short-range structure of each NN interaction. Let us also note 
that the results of Ref. [49] within the two-body case are about 5 — 10 MeV more repulsive 
than ours. 

Once more, we would like to stress the fact that the SCGF approach, complemented 
with the LW formalism, yields thermodynamically consistent results. To this end, we show 
in Fig. 8 the microscopic chemical potential, p,, together with the macroscopic one, /i, as 
a function of density for two temperatures, T = 10 and 20 MeV. Left panels correspond 
to SCGF results for both the CDBONN potential (upper panel) and the Argonne V18 
interaction (lower panel). For the two temperatures, there is a good agreement between 
the microscopic chemical potential, obtained from the normalization condition of Eq. (8), 
and the macroscopic chemical potential, coming from the numerical derivative of the free 
energy density, Eq. (9). For the latter, a centered two-point formula has been used. Let 
us stress that this agreement confirms a posteriori the good behavior of the dynamical 
quasi-particle entropy as a function of the density and also the negligible role of the S' 
term. The approximations involved in the calculation of Sdq do not spoil the consistency 
of the ladder approximation. Results for jl and /i within the BHF approach are presented 
in the right panels. Both chemical potentials agree at low {p < 0.08 fm~^) densities for 
both temperatures. Above this density, however, discrepancies appear due to the increasing 
importance of the rearrangement contribution to the self-energy [51]. At 0.16 fm~^, the 
difference is of ~ 4 — 5 MeV for both potentials and temperatures, and it becomes as large 
as 15 — 20 MeV at 0.30 fm~^. These differences show the lack of consistency of the BHF 
approach at finite temperature, even though the effect is smaller than in nuclear matter [21]. 

The EoS of neutron matter is shown for different temperatures in Fig. 9. For the SCGF 
approach this quantity is computed from the thermodynamical relation p = p (/i — ^) , with 
jl the microscopic chemical potential obtained from Eq. (8). Note that in thermodynamically 
non-consistent approaches the pressure has to be computed from numerical derivatives of the 
free energy with respect to the density. Once again, a remarkable agreement with FP is found 
for the Argonne V18 results, while the CD-BONN interaction leads to a softer EoS. In the low 
density regime, both results agree well with the virial expansion. The effect of temperature 
decreases as density increases and eventually the curves for different temperatures seem to 
collapse to a single (density dependent) value, as expected from degeneracy arguments. This 
high density regime, however, will be mostly affected by the inclusion of three-body forces 
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in the calculations. 

So far, we have plotted all our results as functions of density. To get a more accurate 
insight on the temperature dependence of the different thermodynamical properties of the 
system, we show in Fig. 10 the energy (left panel), entropy (central panel) and free energy 
(right panel) per particle as a function of temperature for a fixed density, p = 0.16 fm~^. 
The results correspond to the Argonne V18 interaction. The agreement between the SCGF 
and the FP energy per particle is confirmed for all temperatures (the FP results have been 
interpolated to this particular density) . The SCGF results are about 2 MeV more repulsive 
than the BHF ones and this repulsive effect is almost temperature-independent. The entropy, 
as expected, is well determined by all the approaches and only some small differences can 
be observed at large temperatures. These small differences, however, are translated into a 
slight disagreements between the FP and SCGF results in the free energy per particle. As 
observed in the right panel, the SCGF results are about 0.5 MeV more attractive than the 
FP ones. The BHF results are about 1 MeV more bound than the SCGF ones. Again, the 
differences between the approaches are rather temperature-independent. This suggests that 
the effect of dynamical correlations on the macroscopic properties are rather insensitive to 
thermal effects. 

It would be interesting to study the effect that sophisticated many-body calculations 
have on the temperature dependence of the different thermodynamical properties. This 
would provide a reliable test for the usually assumed quadratic (linear) temperature de- 
pendences for the energy (entropy). A detailed study of these dependences would however 
need of reliable extrapolations to the low-temperature regime, which in our present ap- 
proach is not possible due to the presence of pairing effects. A thorough analysis of this 
low-temperature regime will be discussed elsewhere. At the moment, using the present data, 
we have parametrized the different thermodynamical quantities in terms of simple fits to 
study the quality of the commonly used approximations. 

The energy per particle of the FFG is well fitted by a quadratic temperature dependence, 
e ~ Co + CgT^, inspired by the Sommerfeld expansion [54]. This expansion is only valid for 
^^1^ << 1, temperatures close to zero. In the fits to the SCGF results, however, we 
have to use the available data between T = 4 and 8 MeV. Fitting a quadratic dependence 
for the energy per particle of the FFG in this temperature range yields a deviation from 
the exact result, Og = — 0.0422 MeV~^, of only 5%. Assuming that this procedure 
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is also valid for the SCGF energies per particle, we find Oe = 0.0339 MeV ^. This 20% 
difference seems too large to be explained simply in terms of the effective mass, which in 
this regime is m* ~ 0.9m. The naive replacement — for instance, does not agree 

with the previous value. The more accurate prediction J;-^2^a "^Im^ t^^]' although closer, is 
also somewhat too large to coincide with the fit to SCGF data. Alternatively, one could 
have tried to obtain an analytic expression for from the low temperature expansion of the 
GMK sum rule formula, but this is difficult due to the non-trivial temperature dependence of 
A{k, cu). Let us also stress that the quadratic thermal dependence of the energy per particle 
comes essentially from the kinetic energy term. The potential energy is rather temperature 
independent and decreases by less than 2 MeV when going from T = 20 to 4 MeV. 

According to Fermi liquid theory, the behavior of the entropy at low temperatures should 
be linear with T. For the FFG, a fit of Eq. (17) in the T = 4 to 8 MeV regime gives a very 
accurate value of = 2ae = 0.0844 MeV~^. A similar one-parameter fit to the SCGF yields 
a slope, Qs ~ 0.0772 MeV~^, in agreement with the Fermi liquid prediction for an effective 
mass, m* ~ 0.92m. This coincides with the value that we obtain for the effective mass at 
k = kp Sit low temperatures. The entropy however shows a clear deviation from this linear 
behavior above 12 MeV. In addition, the FFG prediction = 2ae is partially violated. 
Finally, a quadratic fit to the SCGF data for the free energy per particle, f — fo + o,fT^, 
leads to Qf — —0.0428 MeV This is a somewhat low value, rather close to the FFG 
prediction. In contrast, the FFG relation = — a/ is not well fulfilled. Moreover, the non- 
hnear behavior of the entropy for T > 12 MeV leads to a non-quadratic behavior of F/A 
above this temperature. A consistency check of these low temperature fits is the relation 
Oe — ~ Qf as well as the fact that the zero-temperature extrapolation of E/A, cq = 14.51, 
and of F/A, /q = 14.49 MeV, do coincide. Note, however, that the accuracy of the fits 
depend on the exact position and the number of points considered at low temperatures and 
these are the most sensitive to numerical uncertainties within our approach. Finally, we 
would like to stress that the convergence of our results down to T = 4 MeV implies that 
Tc < 4 MeV. 
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V. CONCLUSIONS 



We have presented the first systematic study of hot neutron matter within the Self- 
Consistent Green's Function formahsm in the ladder approximation for two realistic NN 
interactions, the CDBONN and the Argonne V18 potentials. The calculations cover a wide 
range of densities and temperatures and show the adequateness of this method to account 
for correlations in the microscopic properties of dense, hot hadronic matter. The effect of 
short range correlations in the thermodynamical properties is correctly described by the 
Luttinger-Ward formalism. 

At the microscopic level, short-range effects are particularly important on the spectral 
functions and are manifested in its low and high energy tails. Our results indicate that 
both the location of the quasi-particle peak and the amount of strength in the energy tails 
change substantially with density. On the contrary, thermal effects are very small and only 
affect the region around the chemical potential. The momentum dependence of the real part 
of the on-shell self-energy in the SCGF approach has been compared to the single-particle 
spectrum in BHF-like descriptions. The propagation of hh pairs in the intermediate states 
has a repulsive effect with respect to BHF results. This difference grows with density and 
is larger for momenta below the Fermi momentum, with maximum differences of ~ 25 MeV 
in the range of densities explored here. In the hole-momentum region, in addition, BHF 
and SCGF show different thermal behaviors, with the latter becoming more attractive as 
temperature increases. 

A careful study of the momentum distribution of the system has also been performed and 
the important interplay between thermal and dynamical correlations has been highlighted. 
These effects are well exemplified by n(0), which is customarily used as a measure of corre- 
lation effects. For a given temperature and decreasing density, the system approaches the 
classical limit and the depletion of the momentum distribution increases. For larger densi- 
ties, closer to the degenerate limit, dynamical correlations play a more important role and 
n{0) decreases with increasing density. In general, correlation effects, as measured by the 
depletion, are larger for Argonne V18 than for CDBONN. 

In general, the SCGF energy per particle is more repulsive than the BHF one, indepen- 
dently of the interaction. The magnitude of these differences is governed by the density and 
by the particular structure of the NN interaction, and it is at most of 5 MeV for the range 
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of densities explored (up to 0.32 fm~^). The sensitivity to the NN interaction within the 
SCGF approach appears to be larger than this, with differences of up to 6 MeV in the same 
range. In contrast, BHF results are relatively potential- independent. In any case, in the low 
density regime the energies for both approaches compare very well with the virial expan- 
sion. In addition, there is a very good agreement between the SCGF results for Argonne 
V18 and those of FP for Urbana V14 above a density of 0.08 fm~^. This is possibly due 
to a cancellation between the potential and the many-body dependence of the energy per 
particle in this regime. 

The entropy has been computed within the dynamical quasi-particle approximation, 
which takes into account the effects of correlations in the width of the quasi-particle peak. 
The discrepancies between different approximations to the entropy are rather small, thus 
revealing that the entropy is not affected by correlations. In general, all the approaches lead 
to somewhat lower values than those predicted by the FFG. The free energy for Argonne VIS 
and CDBONN shows substantial differences at large densities due to the different structures 
of the potential. The free energy obtained from the GMK energy and the DQ entropy leads 
to a thermodynamical consistent result, with a good agreement between the microscopic 
and the macroscopic chemical potentials. The differences for the BHF approach can be as 
large as 20 MeV, although in general they are less important than for the nuclear matter 
case. The EoS, which has been computed in a wide range of densities and temperatures, also 
shows a similar potential dependence. In the low density regime, however, all the thermo- 
dynamical quantities show a very good agreement with the virial expansion. The stability 
of our results in this regime shows the robustness of the numerical techniques involved in 
the calculations. 

The temperature dependence of the energy, the entropy and the free energy has also been 
explored. A quadratic dependence is compatible for the energy per particle at p = 0.16 fm~^ 
for the SCGF and BHF approaches. The entropy is only proportional to the temperature 
below T ~ 10 MeV, which in turn translates into a non-quadratic temperature dependence of 
the free energy above this temperature. Moreover, the differences in energy and free energy 
between BHF and SCGF remain constant with temperature, indicating that the effect of 
temperature on dynamical correlations is rather small. In addition, the convergence of the 
results down to T = 4 MeV indicate that no superfluidity appears above this temperature. 
In conclusion, the calculations performed show the potential of the SCGF method to describe 
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accurately the properties of dense and hot matter. The inclusion of pairing effects and three- 
body forces within this formalism will improve the predictions for the microscopic and the 
bulk properties and will provide a very complete description of neutron star matter. 
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FIG. 1: (Color online) Density (left panels) and temperature (right panels) dependence of the 
spectral function, A{k,Lo\ as a function of energy for three different momenta: = (upper), 
k = kp (central) and k = 2kF (lower), with kp the Fermi momenta associated to each density. 
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FIG. 2: (Color online) Real part of the on-shell self-energy for the SCGF and BHF approximations 
at T = 5 MeV (solid lines) and T = 20 MeV (dashed lines). The three panels correspond to densities 
0.08 fm-3 (left), 0.16 fm^^ (central) and 0.24 frn'^ (right). The arrows show the associated Fermi 
momenta. 
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FIG. 3: (Color online) Density (top panels) and temperature (lower panels) dependence of the 
momentum distribution. The results are obtained with the CDBONN (left panels) and Argonne 
V18 (right panels) interactions. 
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FIG. 4: (Color online) Occupation of the lowest momentum state as a function of density for four 
different temperatures. CDBONN (Argonne V18) results are shown on the left (right) panel. The 
double-dotted-dashed line corresponds to the Free Fermi Gas at T = 5 MeV. 
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FIG. 5: (Color online) Energy per particle as a function of density for different many-body ap- 
proaches: SCGF (circles), BHF (triangles), variational (crosses) and virial (solid). The results 
are presented at T = 10 MeV and T = 20 MeV. SCGF and BHF results are computed with the 
CDBONN (left panel) and Argonne V18 (right panel) interactions. The inset in the right panel 
spans the low-density regime for Argonne V18. 
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FIG. 6: (Color online) Entropy per particle as a function of density for different approximations: 
SCGF (circles), BHF (triangles), FP (crosses), Sa (dash-dotted), FFG (double-dot dashed) and 
virial (sohd). The results are presented at T = 10 MeV and T = 20 MeV. SCGF and BHF results 
are computed with the CDBONN (left panel) and Argonne V18 (right panel) interactions. The 
inset in the right panel spans the low-density regime for Argonne V18. 
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FIG. 7: (Color online) Free energy per particle as a function of density for different temperatures 
and different approximations: SCGF (circles), FP (crosses) and virial (solid). SCGF results are 
computed with the CDBONN (left panel) and Argonne V18 (right panel) interactions. The inset 
in the right panel spans the low-density regime for Argonne V18. 



31 



SCGF 



BHF 



> 



• 1— I 
■(-» 

a 

-I— > 
O 
Oh 

13 
O 

• 1—1 

B 

U 



> 



• 1— I 
■(-> 

■I— I 
O 
Oh 

o 

• ^ 

B 
U 



60 

40 

20 



-20 

-40 
60 

40 

20 



-20 

-40 



' 1 ' 1 

- CDBONN 


' 1 1 








• ILi _ 


• J* 


— ILi - 


7t=20 MeV 
? . 1 . 1 


. 1 , 



' 1 ' 1 

- Argonne VI 8 


1 1 jT 


_T=10 MqYj^^ 






— _ 

. 1 , 


-/T=20 MeV 
. 1 . 1 



' 1 ' 1 
- CDBONN 

-* ^ ' V ^ -1-1-1-1 


' 1 ' 










-T=10MeV^ 










A |Ll _ 




— - 


7t=20 MeV 
f . 1 , 1 


, 1 . 



60 
40 
20 


-20 



> 



' 1— I 
■(-» 

a 

■(— > 
O 
Oh 

o 



0.08 0.16 0.24 
P [fm"'] 




0.32 0.08 0.16 0.24 0.32 
P [fm"'] 



FIG. 8: (Color online) Microscopic (symbols) and macroscopic (solid lines) chemical potentials as 
a function of density for temperatures T = 10 and 20 MeV. The left (right) panels correspond to 
the SCGF (BHF) results, while the upper (lower) panels correspond to the CDBONN (Argonne 
V18) interactions. 
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FIG. 9: (Color online) Pressure as a function of density for different temperatures and different 
approximations: SCGF (circles), FP (crosses) and virial (solid). SCGF results are computed with 
the CDBONN (left panel) and Argonne V18 (right panel) interactions. The inset in the right panel 
spans the low-density regime for Argonne V18. 
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FIG. 10: (Color online) Thermodynamical properties as a function of temperature for a fixed 
density of p = 0.16 fm^'^. The SCGF (circles), BHF (triangles) and variational (crosses) results 
are displayed. The dotted lines correspond to fits based on the Sommerfeld expansion. 
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